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C.  I 


Suinmary 

In  their  paper  [1],  Ansarl  and  Bradley  discussed  a 
two-sample  rank  test  for  dispersions  and  suggested  the 
desirability  of  extending  their  results  to  the  prohlem  of 
several  samples.   In  this  paper,  besides  generalizing  their 
results,  we  provide  a  few  additional  non-parametric  tests, 
which  include,  among  others,  the  multi-sample  analogues  of 
the  two-sample  normal  scores  test  of  dispersion  and  the  tests 
considered  by  Mood  [6],  Klotz  [4],  and  Slegel  and  Tukey  [10]. 
The  asymptotic  distributions  of  the  proposed  test  statistics 
are  derived  by  an  application  of  the  author's  theorem  [?]. 
The  asymptotic  efficiencies  of  these  tests  relative  to  one 
another  and  the  ^   test  ([9]  PP-  83-87)  are  computed  in  the 
standard  fashion  along  the  lines  of  the  author's  paper  [7]' 


1.   Introduction 

Let  X.  .(  j=l, . . . ,m. ;1=1, .  .  .  ,c)   be  Independent  samples 

from  populations  with  continuous  cumulative  distribution 

function  F.(x)  =  F(0.(x-v),   where  0.      and   v   are  real 

numbers.  9.      Is  the  scale  parameter  for  F.   and   v,   the 

common  median  for  all   F. ,   which  without  loss  of  generality 

can  be  taken  to  be  zero.   We  are  Interested  In  testing  the 

hypothesis   H  :    9^    =    .  .  .    =    9        against  H-,  :    9 .    ^   9  .      for 
•^  ^  ol  c^         ll'^J 

some  pair   ( 1, j  )  . 

Let  Z^  •  =  Ij   if  the  1th  smallest  observation  from 
the  combined  sample  of  size   N  =  '^   m.   Is  from  the  jth  sample 
and  otherwise  let   Z^  .  =  0.   Then  we  propose  to  consider  the 
following  test  statistics. 


A, 


£(B)  =48^„.(I^_j  -B^^/ 


where 


(1.1) 
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where 


(1.3) 
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and  where  V    < 


.(N) 


<  V^""'   is  an  ordered  sample  of  size 
N   from  a  distribution  Y   and  E   denotes  the  expectation. 
( /5t,t  ■  ,U^J.^    -jCt,.  .   and   A,,   are  normalizing  constants  to  "be 
defined  helow) . 

We  may  note  that  with   c  =  2,   the  statistics   =C(B), 
£(M)   and  r.(J)   where  ^  is  the  standard  normal  distribution 
function,  reduce  respectively  to  the  two-sided  Ansari-Bradley- 
Freund  W  statistic  [1],  Mood's  M  statistic  [6]  and  the  normal 
scores  statistic  [^J?  ^^"^   testing  the  equality  of  dispersions 
of  the  two  populations. 

Reference  to  prior  work  on  the  two-sample  non-parametric 
tests  for  dispersions  may  also  be  found  in  Lehmann  [5],  Sukhatme 
([11],  [12]),  Barton  and  David  [2],  Siegel  and  Tukey  [10], 
and  Klotz  [4],  among  others.   Ansari  and  Bradley  [l]  have 
shown  that  their  W  test  is  equivalent  to  one  independently 
proposed  by  Barton  and  David  [2],   Klotz  [4]  has  shown  that 
the  W  test  is  equivalent  to  one  proposed  by  Siegel  and  Tukey 
[10].   Thus  the  £(B)  test,  defined  above,  may  also  be  regarded 


as  a  generalization  of  the  two-sided  Ansari-Bradley-Barton- 
David-Siegel-Tukey  test  of  dispersion.   Lehmann's  test  is  not 
a  distribution-free  test  and  so  the  relative  efficiencies  of 
this  test  are  not  known.   Sukhatme's  second  test  [12]  is  an 
improvement  over  the  Ansari-Bradley  ¥  test  and  its  generalization 
to  several  samples  merit  investigation;  hut  since  this  does 
not  appear  to  he  possible  with  the  methods  of  the  present  paper, 
we  do  not  intend  to  discuss  this  test  here.   We  shall  confine 
our  attention  to  the  statistics  proposed  in  A,  B,  and  C  above. 
(The  author  is  not  aware  of  any  multi-sample  rank  test  for 
dispersions  in  the  literature). 

In  the  passing,  we  may  remark  that  the  statistics 
proposed  above  may  be  put  in  the  general  framework  of  the 
statistics  £  defined  as 


c 
(1.4)  £  =  2Z  ^ 


J 


<^N,J  -^'N.j'/'^N 
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whe 


re   [i...  .   and  A^^   are  normalizing  constants,  and 


N        f  .. 
(  J  ) 


<l-5'  r^fu.^'l^^.l^nA 


With  Ej^  1  =  I  +  -^  -  1-|  +  ^  -  ^1  '   the  statistic   £ 


reduces  to  the  £(B)  statistic.   With  E^  .  =  (  ^  -  2N~  /    ^^ 
obtain  the  £(M)  statistic  and  with  E^^  ^  =  E^  i  =  -^^fK   J 
we  obtain  the  £(Y)  statistic.   The  tests  "based  on  the  statistics 
£(B),  J:(M)  and  £("?)  will  he  referred  as  the  £  tests. 

The  statistics  Z,    defined  by  (1.4)  with  weights 
E,.  .   different  than  the  ones  proposed  above,  were  suggested 
by  the  author  [7]  for  testing  the  equality  of  location  parameters 
of  the   c  probability  distributions.   Here  we  shall  find  the 
limiting  distributions  of  the  statistics  £  (defined  by  A,  B, 
and  C)  when  the  distributions   F.(x)   differ  only  by  scale 
parameters.   The  formulation  of  the  problem  as  given  in  this 
paper  is  the  same  as  given  in  the  earlier  paper.   However, 
the  earlier  paper  treats  the  problem  somewhat  more  generally, 
and  the  reader  is  referred  to  it  for  additional  background 
and  motivation. 

2.   The  limiting  distributions  of  £  under  local  alternatives. 

To  obtain  the  large  sample  distribution  of  the  statistics 
£,  we  begin  by  determining  the  asymptotic  distributions  of  the 
statistics  T^  .   defined  by  (I.5).   This  is  given  by  the 
following  theorem,  where  the  sample  sizes  m.   are  assumed  to 
tend  to  infinity  in  such  a  way  that  m.  =  s.-n,  n  ->  00;  i=l,...,c, 

Theorem  2.1.   For  n  =  1,2,...   let  X.  .  ( j=l , . . . ,m. ( n ) ; 

i=l,...,c)   be  independently  distributed  according  to 

( n) 
F ( 0 .   X ) .   Suppose  that  the  sequence  of  parameter  points 


)(^)   =    {9['^K...,fMh      satisfy 


(2.1)     -  e[^^  ==  1  +  v./yn 


Let  the  assumptions  of  theorem  6.2  _of  [7]  be  satisfied. 
Then  the  variables   (W-,  ,  ,  .  .  ,W  )   given  by 


have  a  joint  asymptotic  normal  distribution  with  zero  means 
and  covarlance  matrix  whose   ( j , j ' ) th   term  Is 

(2.5)  (  .jj,  --g  JA^ 


where  the   5...   are  the  kronecker  deltas,  and,  where 
J  ,J ' 


(2.4)    A^  = 


i       (Y  ^(x))^dx  -j[       (^"^(x))^dx  J    for  the  £(^)  test. 


Jh   for  the  j:(B)  test, 


•^^Q  for  the  £(M)  test. 


Note  that,  under  the  assumptions  of  Theorem  2.1, 


we  have  (cf .  [7]  ) 


(2.5)  ^^„,j(«'"')  - 


/.:  ['■' 


(H(x)) 


dF.(x)  for  the  Z{f)    test. 


where 


o 


H(x)dP.(x)  +  /    (l-H(x)  )dF  (x) 


for  the  £(B)  test. 


+00 


(r{x)    -  ^)      dF.(x)  for  the  £(M)  test, 


(2.6) 


F  Ax)    =   f(x(1  +  V  //e) 


and 
(2.7) 


H(x)  =  2Z  s,F.(x)  /^ 
1=1   ^  ^    /  a=l 


a 


The  proof  of  the  theorem  follows  as  did  theorem  7.1 
in  [7]  from  theorem  6.2  of  that  paper.   Now  making  the  analysis 
of  variance  transformation 

So  =  ^  ej/2  W.^/a  where   e.  ,  =  ^±/^  H 
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5.  =  y         a.  .  ,¥. 

1   fr^  11'  1 


/a    ;    1=1, . . . , c-1 


where  the   a's   are  chosen  to  make  the  transformation  orthogonal 
and  proceeding  as  In  [j],    we  arrive  at  the  following: 

Theorem  2.2.   Suppose  that  the  assumptions  of  Theorem  2.1 
are  satisfied. 
(a)   If 


(2.i 


11m    v/mT  (  C,,    .  ( 
n->co         J    V   N,j 


,(n) 


-    Cn,j(1W   /a 


exists  and  Is  finite,  then,  for  n  ->  oo,   the  statistic  cC(Y) 
has  a  limiting  non-central  chl-square  distribution  with   c-1 
degrees  of  freedom  and  non-centrallty  parameter 


(2.9) 


1' 


^ 


J-1 


1  Im  ym . 
n->oo 


=N,j('-'""') 


S,J<1' 


/A' 


(b)   If 


(2.10)     li-^^(-«M.j(«'"''  -^N,j(l') 


exists  and  Is  finite,  then,  for  n  ->  oo,   the  statistic   cC(B) 


has  a  limiting  non-centrallty  chl-square  distribution  with 


c-1  degrees  of  freedom  and  non-centrallty  parameter 


(2.11)     A^  =  48  ;^ 


B 


J=l 


u.^y^(^,^.(.>',  .5,,,(i,) 


2 


(c)   If 
(2.12) 


lim  ynT 


(A 


n) 


-  ^L         (1 


exists  and  Is  finite,  then,  for  n  ->  oo,   the  statistic  £(M) 
has  a  limiting  non-central  chl-square  distrlloution  with  c-1 
degrees  of  freedom  and  non-centrallty  parameter 


(2.15)    Aj^  =  180  > 


j=l 


"l^(^.J<'^'"''-*N,j(l') 


Corollary.   If,  in  addition  to  the  assumptions  of 

Theorem  2.2,  the  regularity  conditions  of  lemma  7 . 2  _of  [ 9 ] 
are  satisfied,  then, 

2  \^ 


+00 


(2.14)   A^  =  2_  s  (v  -V) 


X  ^<^~^[F(x)]>  dF(x; 


(2.15)  ^T,  =  ^^Yl   ^j^^j"^^"/  /    ^  f''(x)dx  -  /   X  f^(x)dx  j 
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r 


(2.16)   A^  =  180  ;^  s,(v  -v)^('  /     X  f(x)[2F(x)-l]dF(x) 

j  =  l   ^   J     W-00 


where   f   Is  the  density  of  F   and 


(2.17) 


V  = 


s.  V.  A"  s 
1  \.'  '^— 


1 


Note:   To  get  A-^,   we  need  a  trivial  modification  of  the 
part  (i)  of  lemma  7.2  of  [9]. 


Proof.   To  get  (2.l4)  from  (2.9),  we  note  that 


J^  Gn,j^^'"'^    -    ^ 


,(n) 


N 


,j'i') 


'n/    /    J 


XI  A  f(x(1+v  /y?T)) 
1=1      •  ^       ^ 


dF^(l+v   //n))    -    Tj    [F(x)]dF(x) 


V  J 


A    (x)B    (x)dF(x) 
n  n 


where 


Y~  A.F(x(i+v./y?r)(i+v  /a)  ) 


c 

1=1 


A    (x)    =  - 
n^  c 


J[F(x)] 


jl  A.F(x(i+v./yH3(i+v^./yK)  J  -  F(x) 


B^(x)    =yK  XI  A^   F(x(l+v^/y?T)(l+Vj//Kj       \  -    F(x) 
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Under  the  assumed  regularity  conditions 

u   =  F(x) 


11m      A    (x)    =  A.  J(u) 

n  du 

n->oo 


11m     B    (x)    =    ( V    -    V  .)    X   f (x' 

^        n  J 

n->oo  " 


where  "v   Is  given  by  (2.17).   Since  differentiation  under 

the  Integral  Is  permitted,  the  proof  follows. 

The  proofs  of  (2.15)  and  (2.l6)  are  exactly  analogous, 
In  particular  when  ^  =  ^,   where  ^  Is  the  standard 

normal  distribution  function, 


2 


(2.18)   A  ^  2  X:  s.(v-v)^/^r^  r'[FU)]f(x)dF(x) 


where   4>   Is  the  density  of  ^. 

3.   Parametric  Case 

In  the  parametric  theory,  a  commonly  used  test  Is 
Bartlett's  H  test  for  the  homogeneity  of  variances.   This  test 
Is  shown  by  Box  to  be  very  sensitive  to  departures  from  normality 
assumptions.   We  give  below  a  very  brief  description  of  an 
approximate  test  based  on  the  analysis  of  variance  of  the 
logarithms  of  the  sample  variances.   For  details,  the  reader 
Is  referred  to  Scheffe  [[9],    pp.  83-87),  from  where  the  following 


12 


portion  Is  taken  with  slight  modification. 

Divide  the  1th  sample  X. .,,...  ,X.m..   Into  J.   sut- 

ll         IX  1 

2 
samples.   Let  m. .   denote  the  size  and   s. .,   the  sample 

variance  of  the  jth  sub-sample  of  the  1th  sample.   Denote 

2 
Y.  .  =  log  s.  .   and  assume  that  y^,      the  kurtosls  measure 

has  the  same  value  for  all  populations.   Then  from  Scheffe 

(cited  above),  the  test  for  the  hypothesis  of  the  equality 

of  variances  Is  based  on  the  statistic 

-^2 


(3.1)         7 


c-i 


where 


—  ;^  y~  V.  .(Y.  .  -  f|.) 
^e   1  V  '-'      '-  ^ 


J. 
V.  .  =  m.  .  -  1,  V.  =  ^        V.  .,  r\.    =  y        V.  .Y.  ./v., 


•  _L_L  a^  t^«^ 


A  _ 

1      ^  1 


denominator  of  '~f     converges  to  2+j^     in  probability  as 
n  ->  00.   (Note  that  we  are  assuming  m.  .  ->  <»  as  n  ->  oo). 
Hence  an  asymptotically  equivalent  statistic  is 

(3.2)         (c-l)t=  ^  v"(^.  -  ^)V(2  +  72) 

i 

which  has  a  limiting  non-central  chi-square  distribution  with 
a  non-centrality  parameter  (after  omitting  computations) 
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(30)       A     '■  ^   '  •   -^^ 


^  =  4z:3.(v.-v)y(2+72) 


4.   Asymptotic  Relative  Efficiency. 

It  Is  well-known  ([?])  that  In  the  situations  we  are 
considering,  the  asymptotic  efficiency  of  one  statistic  relative 
to  another  Is  equal  to  the  ratio  of  their  non-centrallty 
parameters.   Hence,  we  have  the  efficiencies  of  the  normal 
scores  £(^)  statistic  relative  to  £(B),  £(M)  and  £("/)  statistics 
as  follows: 


^^•2)  ^£($),£(M)  -  V^M 


(40)  ^£(i),f  =  yv 

It  may  be  remarked  that  the  above  results  agree  with 
the  results  obtained  by  Ansari  and  Bradley  [1],  and  Klotz  [4] 
for  the  two- sample  case,  and  hence  the  results  of  this  paper 
as  well  as  of  these  authors  apply  directly  to  the  c-sample 
problem.   We  would,  in  particular,  draw  the  attention  of  the 
reader  to  Section  3  of  [h]    and  Section  7  of  [l]  where 


the  efficiency  comparisons  for  different  densities  of  the 
£(J),  £.(M),  £(B)  and  4  tests  are  made  for  the  two-sample 
problem  and  which  are  shown  here  to  he  valid  also  for  the 
c- sample  problem. 


Ih 
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